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Ca, the lightest double beta decay candidate, is the only 
one simple enough to be treated exactly in the nuclear shell 
model. Thus, the /3j3(2v) half-life measurement, reported 
here, provides a unique test of the nuclear physics involved 
in the (3(3 matrix element calculation. Enriched 48 Ca sources 
of two different thicknesses have been exposed in a time pro- 
jection chamber, and yield Tff 2 = (4.3tj' 4 [stat.]±1.4[syst.]) x 
10 years, compatible with the shell model calculations. 

PACS numbers: 21.10.Tg, 21.60.Cs, 23.40.Hc, 27.40. +z, 
23.40.-s 

Among all candidates for double beta decay, 
48 Ca— > 48 Ti is unique, since it is the only one which can 
be treated "exactly" in the nuclear shell model by solv- 
ing the problem of eight nucleons distributed within the 
fp shell without truncation. Consequently, this decay 
has been a favored testing ground of nuclear theories [Q . 
However, until now only a lower limit of the 2v decay 
half-life, T 2 ^ 2 > 3.6 x 10 19 years, has been determined 
experimentally gj. 

In a recent paper || it was shown that the nuclear shell 
model, constrained by the requirement that it describes 
well the spectroscopy of the A = 48 nuclei, restricts the 
corresponding half-life also from above, T 2 ^ 2 < 10 20 y. 
Therefore the experimental observation of the decay ac- 
quired added significance; if it turns out that the shell 
model cannot predict this theoretically tractable rate, 
we have to wonder about our ability to describe the nu- 
clear matrix elements in the more complex nuclei, includ- 
ing perhaps even the matrix elements responsible for the 
neutrinoless decay. Here we will present the result from a 
new experiment using a time projection chamber (TPC). 

The search for double beta decay of 48 Ca was among 
the first to be attempted in live-time experiments be- 
ginning in the early fifties (for an extensive chronology 
see With the largest energy release of all (3(3 candi- 
dates 48 Ca = 4.3 MeV) has a (3(3{2v) sum-energy 
spectrum that extends to higher energies than most ra- 
dioactive background. Yet calcium has a tendency to 
harbor chemically similar radio-impurities such as 90 Sr 
and 226 Ra, which do intrude on a major fraction of the 
48 Ca spectral range. 

When the two (3 particles are tracked in a TPC they 
are seen to both carry negative charge, originate from 
a common point, and have separately measured ener- 



gies. Although this distinctive visualization eliminates 
the bulk of unrelated activity, there remain several well- 
known mechanisms for production of negative electron 
pairs that constitute background for (3(3 decay. The most 
serious of these are Moller scattering of single (3 particles, 
and cascades in which a y-ray internally converts or 
Compton scatters. These processes are fed principally 
by the primordial decay chains, but also by cosmogenic 
and man-made radionuclides. Decay-chain induced back- 
ground events can often be tagged by a or (3 particles 
from neighboring links if the source is thin enough to 
allow the tagging particles escape into the TPC gas. 

The miniscule 0.187% natural abundance of 48 Ca 
makes enrichment both necessary and expensive. Po- 
tential exposure to loss of costly isotope is a deterrent 
to chemical purification or conversion to the lightest sta- 
ble compound. The material used was supplied by the 
Kurchatov Institute as finely powdered CaC03 enriched 
to 73% in 48 Ca, and relatively free of U and Th (< 0.8 
ppb by mass spectroscopic analysis). The speculative 
assumption of secular equilibrium and the absence of se- 
rious amounts of other high Q- value impurities, led to a 
projected experimental (3(3 sensitivity in excess of 10 20 
years for this material. Accordingly, a relatively thick 
(3(3 source was made to maximize the number of exposed 
48 Ca nuclei. The CaC03 powder was injected into a large 
glass box in bursts of compressed gas, allowed to settle 
uniformly onto a 4 fj,m Mylar substrate, then fixed with a 
mist of Formvar. Two such deposits, face-to-face, formed 
the first (3(3 source, with a total of 42.2 g of CaC0 3 (18.5 
mg/cm 2 total thickness with substrate and binder.) 

The shielded UC Irvine TPC |,§ containing the (3(3 
source as the central electrode in a magnetic field was lo- 
cated in a tunnel at the Hoover Dam under a minimum of 
72 m of rock. Data were recorded on magnetic tape, and 
subsequently passed through stripping software to select 
clean le~ and 2e~ events. The 2e~ events were individ- 
ually scanned by a physicist. All unambiguous negative 
pairs emitted from opposite sides of the source with a 
common point of origin were fitted with helices, and the 
parameters written to a (3(3 candidate file or a 214 Bi file, 
depending on whether a 214 Po a particle appeared at the 
vertex within the following millisecond. The far more nu- 
merous le _ events were fitted automatically by software 
and also written to a parameter file. 

The lone electron (le~) spectrum plotted against ki- 
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netic energy (K) in Fig. [j] represents the total beta ac- 
tivity of source contaminants, and can be broken down 
into the contributions from individual radionuclides by a 
least squares fit. 90 Sr (2250 /iBq/g), 226 Ra (530 /iBq/g) 
and their daughters account for the bulk of the spectrum. 
Contributions from 137 Cs (940 /iBq/g), and daughters of 
228 Ra (90 /iBq/g), are also present. The two Ra activi- 
ties being much larger than the mass spectroscopic limits 
on U and Th, indicate severe breaking of equilibrium in 
the respective series. The daughters of greatest concern 
(Qp > 2 MeV) are 90 Y, 214 Bi, 228 Ac, 212 Bi, and 208 T1. As 
a check of the fitting procedure, the energy spectrum of 
electrons tagged by 214 Po a-particles was noted to closely 
match the fitted 214 Bi component when adjusted for the 
a escape probability (P Q = 0.24 ± 0.01 from an indepen- 
dent measurement.) 
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FIG. 1. The lone electron spectrum for the thick source, 
and the most important of the fitted components. From top 
to bottom at low energy are the total spectrum, 90 Sr, 226 Ra, 
and 228 Ra. Daughters are assumed in equilibrium, with the 
exception of 210 Pb. Its daughter 210 Bi was fitted separately, 
then included in the 226 Ra curve. 

The activities determined from the lone electron spec- 
trum were used as input to a Monte Carlo calculation 
of the 2e~ background. 90 Sr and its daughter 90 Y are 
essentially pure (3 emitters, and contribute only through 
Moller scattering. The other high nuclei have com- 
plex decay schemes with multiple gamma rays [Q, all of 
which were included in the Monte Carlo with their cor- 
responding conversion coefficients. Most of these sim- 
ulated 2e~ backgrounds were directly testable against 
TPC measurements, with good agreement: The mea- 
sured 21 Bi component was simply the 2e _ data sub- 
set tagged by the 214 Po a, and corrected for the a es- 
cape probability. The 90 Y 2e~ measurement was pro- 
vided by a drop of 90 Sr solution applied to a natural iso- 
topic replica of the 48 Ca source, and placed in the TPC. 
The 212 Bi and 208 T1 2e~ measurements were scaled from 
those produced by an injection of 220 Rn, by comparing 



observed rates of the rapid 212 Bi- 212 Po, (3-a sequence. 

Since the Monte Carlo 2e~ rates were derived from 
intrinsic activity levels in /iBq/g, their agreement with 
direct measurements also confirms the Monte-Carlo pre- 
dicted 2e~ efficiency of the TPC. A Monte-Carlo gener- 
ated background spectrum was essential only for 228 Ac 
where we have no TPC measurement, but in our 2e~ 
background model we elected to use the smoother, better- 
statistics Monte Carlo spectra for the other contributions 
as well. 

An alternative determination of 2e~ background was 
carried out by a separate subgroup of the collaboration, 
and included independent Monte Carlo calculations and 
a greater reliance on the above-mentioned TPC measure- 
ments as opposed to lone-electron fits. We refer to this 
direct measurement method as analysis "A", and the 
lone-electron based method as analysis "B". To elimi- 
nate events with the poorest energy resolution, analysis 
B included a cut on electrons making the smallest angles 
with the magnetic field, (|cos(0) | < 0.9). 

The 2e~ sum spectrum from the 48 Ca source, after 
event-by-event removal of the a-tagged 214 Bi, is shown in 
Fig. ||a with a singles threshold of 400 keV. The various 
remaining background spectra as determined by analy- 
sis A are superimposed. The residual spectrum following 
background subtraction, appears in Fig. ^o. 
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FIG. 2. (a) The measured 2e~ sum-energy spectrum and 
background spectra for the thick source with a 400 keV singles 
threshold. From top to bottom at low energy are the total 
measured spectrum, Moller events, 214 Bi, and daughters of 
228 Ra. (b) Residual /3/3 candidates and Monte Carlo 0/3 spec- 



trum at the corresponding T-wj 
this threshold (dashed line). 



4.1 x 10 y obtained at 



The above procedure was repeated for a series of singles 
thresholds and one relatively high sum threshold, with 
consistent results. The corresponding half-lives were cal- 
culated in each case, as shown for two of the threshold 
combinations in Table Q. The 2 MeV sum threshold was 
accompanied by an additional singles spectrum cut on 
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the strong 214 Bi conversion line at 1.3 MeV. Quoted er- 
rors in the table are statistical. The calculated half-life 
is independent of threshold and analysis A or B within 
errors. 

Since the half-life precision was degraded by the unex- 
pectedly large Ra and Sr contamination, the 42.2 g source 
was replaced after 2440 hours exposure, with one contain- 
ing only 10.3 g of enriched CaCC>3, and total thickness 
5.4 mg/cm 2 . The thinner source was exposed 4001 hours, 
and these results are also included in Table |. The 214 Bi 
rejection improved greatly as a result of an increase in Pq, 
from 0.24 to 0.69, and the 90 Y component was weakened 
by the large reduction in target mass for Moller scatter- 
ing. 

In view of these significant improvements, the rela- 
tively small increase in signal-to-background ratio ob- 
served for the thin source was disappointing. We now 
understand that multiple Coulomb scattering caused a 
nonlinear dependence of total Moller cross section on 
source thickness. In addition, the thick source degraded 
a larger percentage of electrons to sub-threshold energies, 
and thereby suppressed background more effectively than 
the generally more energetic (3(3 events. Nevertheless, 
it is encouraging that despite different dependence on 
source mass for the 48 Ca signal and the background rates 
(particularly 214 Bi,) the half-lives derived from thick and 
thin sources agree within statistics. 

The thick and thin source residual spectra were each 
corrected for energy loss and efficiency distortion, then 
combined for the Kurie plot shown in Fig. ||. Includ- 
ing the singles threshold (e) in an integration of the 
Primakoff-Rosen spectral approximation Q| produced a 
Kurie formula that retains the linear energy dependence 
of a transformed (3(3{2u) spectrum in the presence of a 
singles threshold cut. The resulting formula is then 



dN 
dK 



where 



/((K-2e)(f (K) + UK))) 
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oc (Qp/3 - K) 



fo(K) = K 4 /30 + K 3 m/3 + 4K 2 m 2 /3 + 2Km 3 + m 4 
and 

f e {K) = e(K-e)(K 2 /l5+2Km/3+2m 2 /3+e(K-e)/5) . 

with electron mass (m) . The Kurie plot energy intercept 
at 4.2 ± 0.1 MeV is consistent with the 48 Ca value 
of 4.271 ± 0.004 MeV ||. The small error bars resulting 
from the Kurie transformation, have been omitted in the 
figure. By comparison, the a-tagged 214 Bi 2e~ events 
produce a distinctly bowed Kurie plot, and the plot for 
measured 208 T1 (Q^ = 5.0 MeV) events is grossly nonlin- 
ear. 

Since the thick and thin sources were exposed sepa- 
rately, corresponding pairs of columns in the table can 



be combined as independent measurements. For exam- 
ple, combining thick and thin (3(3 events for 0.400/0.800 
MeV singles/sum thresholds, analysis A, yields T 2 J 2 = 

(4.3t 2 ; 4 ) x 10 19 y. Either of the other two thick-thin pairs 
would combine to give an equally valid result. However, 
since the other two results would not be statistically in- 
dependent from the first, we do not attempt a grand 
average. Rather we choose the above number, since it 
includes a broader range of the spectrum than the higher 
sum threshold, and we include the A-B analysis difference 
in the systematic error. The remainder of the systematic 
error is largely in the detector efficiency. Thus we quote a 
final result of T 2v /2 = (4.3± 2 ; 4 [stat.] ± 1.4[syst.]) x 10 19 y. 
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FIG. 3. High-energy portion of the Kurie plot for residual 
(3(3 candidates combined from thick and thin sources at a sin- 
gles threshold of 400 keV. (F(K,e) = (K-2e)(f (K)+f e (K)). 
See text.) 

This observation is consistent with the spectroscopy- 
constrained shell model of reference [||. For most (3(3 
nuclei the agreement between the experimental half-life 
and the calculated one has been obtained only with the 
empirical quenching of the Gamow- Teller (GT) strength, 
i.e., the multiplication of all matrix elements of the op- 
erator <tt± by a universal quenching factor 0.77. Taking 
into account the range of lifetimes compatible with the 
limits of Ref . Q , and the experimental uncertainties de- 
scribed above, we cannot, alas, make a definitive state- 
ment about the GT quenching for the 2v double beta 
decay (but for the optimal shell model the quenching is 
clearly required). On the other hand, some calculations 
fTpf that are unconstrained by the spectroscopy test, sug- 
gest a much longer half- life for 48 Ca in disagreement with 
the present experiment. 

Although the 48 Ca enriched sample contained traces of 
radium and strontium activity, backgrounds were well de- 
fined by associated alpha particles and the lone electron 
spectrum. The residual data give half-life values that are 
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TABLE I. Breakdown of counts from the two sources for two energy thresholds. 



Thick Source (0.0775 moly) Thin Source (0.0310 mol-y) 



Singles threshold (MeV) 0.400 0.200 0.400 0.200 



Sum threshold (MeV) 
214 Bi 1.3 MeV line cut 


0.800 
No 


2.000 

Yes 


0.800 
No 


2.000 
Yes 




Analysis A 


Analysis B 




Analysis A 


Analysis B 




Small polar angles cut 


No 


Yes 


No 


No 


Yes 


No 


Counts 














With a tag 


55 


50 


6.8 ± 1.4 a 


79 


72 


8.0 ± 1.6 a 


Without a tag 


500 


472 


72 


142 


134 


21 


Backgrounds b 














Untagged 214 Bi 


189.1 ±25.5 


186.3 ±35.7 


23.4 ±4.9 


40.2 ±4.5 


36.4 ±8.9 


4.1 ±0.8 


Moller events c 


191.5 ±4.6 


163.4 ±5.6 


3.2 ±0.6 


53.7 ± 1.7 


46.1 ± 2.1 


1.5 ±0.3 


228 Ac, 212 Bi, 208 T1 d 


42.8 ±7.5 


29.4 ±5.6 


4.0 ±0.8 


19.8 ±3.5 


13.3 ±2.1 


1.4 ±0.3 


Total background 


423.3 ± 27.0 


379.1 ±36.6 


30.7 ±5.0 


113.7 ±6.0 


95.8 ±9.4 


7.0 ± 0.9 


PP events 76.7 ± 35.0 
Efficiency c 0.0973 
T?/ 2 (10 19 y) 4.1±?;1 
Signal/Background 0.18 
Kurie plot intercept (MeV) 4.1 ± 0.1 


92.9 ±42.5 
0.0902 

o 0+2.7 
J - z -l.() 

0.25 


41.3 ±9.8 
0.0414 

o q+1-0 

a -°-0.6 

1.35 


28.3 ± 13.3 
0.1046 
4 8+ 4 ' 2 

0.25 
4.4 ±0.1 


38.2 ± 14.9 
0.0972 

o 0+2.1 
°-°-0.9 

0.40 


14.0 ±4.7 
0.0469 
4.3l 2 ; 2 
2.04 



a Scaled from counts at lower threshold by the ratio observed in a larger 214 Bi data set. 
b See text. 

c Exclusive of 214 Bi which is included in the row above. 
d Exclusive of Moller,which is included in the row above. 
°From Monte Carlo simulation. 



consistent between two (3/3 sources of different thickness, 
and among various energy thresholds. The correspond- 
ing Kurie plots all intercept the energy axis near the 48 Ca 
Q/3/3 value, unlike plots made from measured samples of 
the various backgrounds. We believe these results con- 
stitute strong evidence for 48 Ca double beta decay at a 
half-life supporting the relatively rigid shell model calcu- 
lations for this lightest double beta decay nucleus. 
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